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Abstract 

The method of multiple scales is used to study the wobbling kink of the cp"^ equation. The 
amplitude of the wobbling is shown to decay very slowly, as and hence the wobbler turns 

out to be an extremely long-lived object. 
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I. INTRODUCTION 



Since the early 1960s, the one-dimensional 0^ theory has been among the most heavily 
utilised models of statistical mechanics and condensed-matter physics [1]. At the same time, 
it served as a testing ground for a variety of ideas in topological quantum field theory 
and cosmology [^. The equation of motion for the model reads 



20« - 2<^-- - + 0' = 0, (1) 

and the fundamental role in applications is played by its kink solution, 

(f){x, t) = tanhx. (2) 

Mathematically, the 0^ kink has a lot in common with its sine-Gordon counterpart; the 
two kinks are the simplest examples of topological solitons in one dimension. There is an 
important difference though; the sine-Gordon equation is integrable whereas the (j)'^ theory 
is not. Accordingly, the kink-antikink interaction becomes a nontrivial matter in the (p^ 
case . Another (not unrelated) difference is that unlike the kink of the sine-Gordon 

equation, the 0^ kink has an internal mode — an extra degree of freedom which allows for 
oscillations in the width of the kink. Although these oscillations are accompanied by the 
emission of radiation (a manifestation of the nonintegrability of the 0^ model), the radiation 
seems to be quite weak and the oscillations are sustained over long periods of time. For small 
oscillation amplitudes, this periodically expanding and contracting kink (termed wobbling 
kink in literature) can be characterised simply as a linear perturbation of the stationary kink 
([2]); however for larger amplitudes one needs a self-consistent fully-nonlinear description. 

The wobbling kink was discovered in the early numerical experiments of Getmanov \^ 
who interpreted it as a bound state of three ordinary kinks. (For a more recent series of 
numerical simulations, see j?!.) Rice and Mele have reobtained this nonlinear excitation 
within a variational approach employing the width of the kink as a dynamical variable 
js]. Segur then constructed the wobbler as a regular perturbation expansion in powers of 
the oscillation amplitude j^. He has calculated the first two orders of the perturbation 
series and noted the likely occurrence of unbounded terms at the third, e^-, order, implying 



the consequent breakdown of the expansion. Subsequently, Sukstanskii and Primak [10 1 



examined the mobility of the wobbler using a variant of the Lindstedt-Poincare technique 
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where the velocity of the kink is adjusted so as to ehminate secular terms at the lowest 
orders of the perturbation expansion. Although this approach provides valuable insights 
at the lowest orders, it is not suitable to deal with the secular terms at the e^-order and 
therefore gives rise to a nonuniform expansion. Finally, Kiselev studied the perturbed 0^ 
kink using the Krylov-Bogoliubov method HI]. His two papers have a wealth of useful 
formulas; however a self-consistent system of equations for the kink's parameters was not 



derived in 



111 ] and the long-term evolution of the wobbling kink has remained unexplored. 



The aim of this paper is to construct the wobbling kink as a singular perturbation expan- 
sion which remains uniform for long times and large distances. Our approach employing a 
sequence of space and time scales will also yield equations for the amplitude of the wobbling 
mode which will be used to analyse its long-term evolution. 



II. MULTISCALE EXPANSION 

Making the change of variables (x, t) where ^ = x — J^v{t')dt', t = t, we 

transform Eq.([T]) to the co-moving frame: 

1 ?7 1 — v'^ 

-(prr - V(p^r - ' (p^ " + 0^ = 0. (3) 



Like the authors of 



lOj, we shall determine the kink's velocity v{t) by imposing the condition 



that the kink be always centred at = [i.e. at a; = v{t')dt']. 

At first glance, the inclusion of the function v{t) is unnecessary: having constructed a 
quiescent wobbling kink, we could make it move simply by a Lorentz boost. By introducing v 
explicitly in Eq.(l3]) we wish to check whether the wobbling kink can drift with a nonconstant 
velocity. The soliton moving with a variable v{t) could obviously not be Lorentz-transformed 
to the rest frame. 

We expand the field about the kink 0o = tanh^: 

= 00 + 601 + 6^02 + . . . . (4) 

Here e is a (formal) small parameter. Substituting (jlj) in ([3]) and setting to zero coefficients 
of like powers of e would constitute Segur's approach which is expected to produce secular 
terms in the expansion ^]. To avoid these, we introduce stretched space and time scales 

X„ = eX T„ = eV, n = 0,l,2,..., (5) 



with the standard short-hand notation for the derivatives: 9„ = Dn = In the 

hmit e ^ 0, the different scales become uncoupled and may be treated as independent 
variables: (pn = (pni^o, ^i, ■■■',To,Ti, ...). We also assume that v is small and slowly varying, 
i.e. V = eV where V = V{Ti,T2, . . .) is of order 1. Expanding ^ and ^ using the chain rule, 
and substituting along with the series (jl]) in Eq. , we equate coefficients of like powers of 



A. Linear corrections 

At O (e^), we obtain the linearisation of Eq.([T]) about the kink 0o = tanhAg: 

^Dl<p, + C<p, = 0, (6) 
where we have introduced the Schrodinger operator 

£ = -^^2 - 1 + = -l^^ + 2 - 3 sech^ Aq. (7) 
The general solution of the variable-coefficient Klein-Gordon equation ([HD can be written as 
01 = CoVoiXo) + Ae*-°^°yi(Ao) + c.c. + 0,j(Ao,To), (8) 

/oo 
[7^(p)e*'^(p)^« + 7^*(-p)e-^-(^')^°] ypiXo)dp. (9) 

Here yo and yi are eigenf unctions of the operator C associated with eigenvalues and |, 
respectively: 

7/o(^o) = sech2Ao; (10) 
l/i(Ao) = sech Aotanh Aq. (11) 

The functions yp(Ao) are solutions pertaining to the continuous spectrum of £: 

Cyp{Xo)= (^2 + ^^yp(Ao); (12) 

these were constructed by Segur j^: 

9.(9. - ir,\ 

(13) 



yp = e^^^« 



^^3(1 — ^tani^^^(i^tg^j^hAo) - 'fV + tanhAo)' 



l+p2 A+p^ 
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The internal mode frequency cuq = while the phonon frequencies uj{p) are given by 
uj{p) = a/4+~p^ > 0. Finally, the coefficients TZ{p) and A are complex; Co is real, and c.c. 
in ([H]) stands for the complex conjugate of the immediately preceding term. 

We set the coefficient of the translation mode to zero as the kink is assumed to be centered 
at Xq = 0. Since we are interested in the wobbling of the kink sustained over long periods of 
time, we also set the radiation amplitudes TZ{p) = 0. As a result, the first-order perturbation 
comprises only the wobbling mode: 

01 = A(Xi, . . . ; Ti, . . .) sech Xq tanh Xoe''^°'^o + c.c. (14) 
B. Quadratic corrections 

At the second order in the perturbation expansion we arrive at the equation 

^Dl<P2 + Ch = F2{Xo,..-To,...), (15) 

where the forcing term is 

F2 = {dodi - DoDi)0i - 3000? + ^^001 + ^I^i^9o0o - ^^'9o'0o. (16a) 
Substituting for 0o and 0i, this becomes 

Fa = -6|Ap sech^ Xq tanh^ Xq + ^DiV sech^ Xq - sech^ Xq tanhXo 
+ [{diA + icuoVA) (2 sech^ Xq - sech Xq) - icuoDiA sech Xq tanh Xq] e^"^"^" + c.c. 

-3A^ sech^ Xo tanh^ XoC^^"^"^" + c.c. (16b) 

The To-independent term in Eq. (116bp and the term proportional to 6*"^°^° are resonant with 
the two discrete eigenmodes of the operator in the left-hand side of f fT^ . while the term 
proportional to e^*"^""^" is resonant with its continuous spectrum. The latter part of the 
forcing is localised in the region near the origin and acts as a source of radiation which 
spreads outward from there. 

Once transients have died out, the e^-correction will consist only of the harmonics present 
in the forcing, i.e. 

02 = ^ + + <^(2)e2^"«^« + C.C, (17) 
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where (^2 \ ^^2^ ^^^1 ip^2^ are functions of Xq which satisfy 

C^f = -el^l^sech^Xotanh^Xo + ^DiV sech^ Xq - 1/^ sech^ Xq tanhXo, (18) 
{C - |)(^J^ = {diA + zcJo^A)(2sech^Xo - sechXo) - iwo/^iAsechXotanhXo, (19) 

and 

(£ - 6)(^^^^ = -3^2 sech^ Xq tanh^ Xq. (20) 

These equations admit bounded solutions if and only if their right-hand sides are orthog- 
onal to the corresponding eigenfunctions of the operator C, Eqs.lITOl) and ( fTTj) . For this to 
be the case, we must set DiV = and DiA = 0. The variation of parameters yields then 

.(0) -0\A 



ip^2 = 2|/l|^sech'XotanhXo + - 3\A\' j Xosech'Xo (21) 
and 

= -{diA + iiUoVA)Xo sechXotanhXo. (22) 

Although the function (p'^^ decays to zero as |Xo| oo, the product e(p^2^ becomes 
greater than the first-order correction ?/i(Xo) for each fixed e and sufficiently large |Xo|. 
Consequently, the term e^02 in the expansion (jl]) becomes greater than the previous term, 
6^01, leading to the nonuniformity of the expansion. In order to obtain a uniform expansion, 
we set 

diA + itooVA = 0, (23) 

which gives A = Aoe-^"°^^S with Aq = Ao{X2,Xs, ...^T^^T^, ...). 

We also note the terms proportional to Xq sech^ Xq in Eq. ( !2T]) . These terms do not grow 
bigger than the previous term, 00 = tanhXo, yet they become larger than the difference 
00 — 1 as Xq — oo and 0o — (—1) as Xq — > — oo. If we attempted to construct the asymptotic 
expansion of the function 0—1 at the right infinity or the function 0+1 at the left infinity, the 
terms in question would cause the nonuniformity of these. Since the function Xq sech^ Xq 
is nothing but the derivative of tanh(fcXo) with respect to fc, these terms represent the 
variation of the kink's width. Hence the potential nonuniformity of the expansion can be 
avoided simply by incorporating them in the variable width [see Eq. ( l36l) below] . 

We now turn to the remaining nonhomogeneous equation, Eq. (j20l) . The variation of 
parameters gives 

^f=AV(Xo), (24) 



/(Xo) = l[6 tanhXo sech^ Xq + (3 - tanh^ Xq + iko tanhXo)( J2 - J^)e''''^'> 
8 

+ (3 - tanh^ Xq - iko tanh Xq) Jsc-'^^^^o] . (25) 
Here the wavenumber k^ = ^/8] the function J2 = J2{Xo; k^) is defined by the integral 

e^^sech^idi, (26) 

-00 

and the constant is the asymptotic value of J2{Xq\ k^) as Xq — > 00. 

The solution ([21]) describes right-moving radiation for positive Xq and left-moving radia- 
tion for negative Xq. The function fl24p is bounded but does not decay to zero as |Xo| — > 00. 



C. The amplitude equation 

Collecting terms of order gives the equation 

^DoV3 + /:03 = i^3, (27a) 

where 

1 



= {dodi - DoDi)<f)2 + {dod2 - DoD2)(j)i + -{df - Dl)<j)i - <f)l - 6<f)o(j)ih + VDodo(j)2 



+ VDod,(j)i + VD^do^i + ^D2Vdo4o - l^'^o'/'i- (27b) 

Having evaluated F3 using the known functions 0o; 0i and 02, we again decompose the 
solution into simple harmonics as we did at O (e^). The solvability condition for the zeroth 
harmonic in equation ( ]27bl) gives D2V = 0, which means that V remains constant up to 
times t ~ e~^. The solvability condition for the first harmonic produces 

t'^D2A + C\A\^A-V^A = 0, (28) 
3 

where 



C = 6 / sech^ Xo tanh^ Xq 



^sech^Xo tanh Xo- 3X0 sech2Xo + /(Xo) dXo. (29) 



The imaginary part of C, can be evaluated analytically: 



Im C = , ° , , = 0.04636, (30) 

^ sinh^ (7rA;o/2) ^ ' 



while the real part can only be obtained numerically: 

ReC = -0.8509. (31) 

Denoting eAo = a{t) the unsealed amplitude of the wobbling mode and recalling that 
V = eV, we express the amplitude equation (^Hj) in terms of the original variables: 

iat = -'^ \a\^a + "^v'a + O (|a|^) . (32) 



We are referring to Eq. fl32|) as the "master" amplitude equation. The master equation 
contains solvability conditions at all orders [which arise simply by expanding the derivative 
d/dt in powers of e] but unlike any particular amplitude equation, it is applicable for all 
times. The master equation is the final product of the asymptotic analysis; all the conclusions 
about the behaviour of the wobbler's amplitude shall be made on the basis of this equation. 
The modulus of a is governed by the equation 

^|a|2 = -u;oImC|a|' + C(|a|^) . (33) 

Since Im^ > 0, the amplitude of the wobbling is monotonically decreasing with time: a 
constant emission of radiation damps the wobbler. The decay law is straightforward from 
dMD: 

|a(t)P = \^ = \^ (34) 

' ^ ^' l + cuoImC|a(0)|2t 1 + 0.08030 X |a(0)|2t' ^ ' 

When a(0) is small, the decay becomes appreciable only after long times t ~ |a(0)|^^. The 
decay is slow; for times t > 12.5 x |a(0)|"^, Eq.(l3ill gives \a\ ~ 

We have verified the above decay law in direct numerical simulations of the full partial dif- 
ferential equation ([T]). As the initial conditions, we took 0(a;, 0) = tanha; + 2ao sech x tanhx 
with some real and 0j(x, 0) = 0. After a short initial transient, the solution was seen to 
settle to the curve flMl) with |a(0)| close to Oq, see FigJH 

The equation (l32ll gives us the leading-order contributions to the frequency of the wob- 
bling: 

uj = ujQ [1_ i^;2+ iReC|a|2 + C(|a|^)] , (35) 

with Re C < as in (I3T1) . The | a \ ^-term here is a nonlinear shift from the linear frequency 
ujq = Vs. The f^-term comes from the transverse Doppler effect. We could have obtained 
this term simply by calculating the wobbling frequency in the rest frame and then multiplying 



the result by the relativistic time-dilation factor y/l — v'^ (which becomes 1 — |f ^ for small 
v). 
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FIG. 1: The decay of the free wobbhng kink. Crosses: \a{t)\ as obtained from the direct numerical 
simulations of Eq.([T]) with the initial conditions in the form (j) = tanh x + 2ao sech x tanh x, (pt = 0, 
with qq = 0.05. Solid curve: equation (j34p with |a(0)| = 0.0492. The inset shows the first 250 
periods of oscillation; during this time the amplitude drops by less than 10%. 

III. CONCLUSION 

One result of this project is the uniform asymptotic expansion of the wobbling kink. In 
terms of the original variables, this expansion reads 

= tanh + asechetanhee^^»(*-'^«) + c.c. 

+ 2|a|2sech2etanhe + aV(Oe''"»(*-"^) + c.c. + (|a|=^) . (36) 

Here ^ = x — vt and /(^) is given by Eq.( l25|) . The first term describes the "background", 
stationary, kink with the width modified by the wobbling. [Note that we have incorporated 
two Xq sech^ Xq terms of the sum fl2ip into the width of the kink.] The second term is the 
wobbling mode itself; the third one accounts for further stationary deformation of the kink's 
shape, and the last term represents the second-harmonic radiation from the wobbler. 

Our second result is the amplitude equation ( l32l) and the conclusion of the extreme 
longevity of the wobbling mode which follows from (!32l) . In view of its anomalously long 
lifetime, the wobbling kink can be regarded as one of the fundamental nonlinear excitations 
of the 0^ theory, on par with the nonoscillatory kinks and breathers. 
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